Mult if r act al properties of aperiodic Ising model: 
role of the geometric fluctuations 
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The role of the geometric fluctuations on the multifractal properties of the local magnetization of 
aperiodic ferromagnetic Ising models on hierachical lattices is investigated. The geometric fluctua- 
tions are introduced by generalized Fibonacci sequences. The local magnetization is evaluated via 
an exact recurrent procedure encompassing a real space renormalization group decimation. The 
symmetries of the local magnetization patterns induced by the aperiodic couplings is found to be 
strongly (weakly) different, with respect to the ones of the corresponding homogeneous systems, 
when the geometric fluctuations are relevant (irrelevant) to change the critical properties of the 
system. At the criticality, the measure defined by the local magnetization is found to exhibit a 
non-trivial F(a) spectra being shifted to higher values of a when relevant geometric fluctuations 
are considered. The critical exponents are found to be related with some special points of the F(a) 
function and agree with previous results obtained by the quite distinct transfer matrix approach. 



I. INTRODUCTION 

The introduction of deterministic aperiodicity in the 
couplings and fields has been intensively used to mimic 
of the effect of disorder on physical homogeneous systems. 
This approach, opposed to the more conventional random 
disorder distribution of coupling constants and fields, has 
been first used in the analysis of electronic systems, and 
is now being used to investigate the critical behavior of 
magnetic systems on Euclidean lattices, as well as on 
fractal hierarchical lattices. 

Hierarchical lattices have been largely used as a frame- 
work for analyzing both homogeneous!] and quenched dis- 
ordered magnetic systems!! It is well known that they 
offer only crude approximations for the critical properties 
of the corresponding homogeneous systems on Euclidean 
spaces with the same fractal dimension, but they provide 
an useful tool to investigate, via exactly solvable models, 
the effect of both random and deterministic disorder. For 
instance, the short range Ising spin glass model defined 
on the diamond hierarchical lattice has been analyzed by 
means of a methodology that encompasses an exact real 
space renormalization group decimation and an exact re- 
cursive procedure to calculate the local magnetization of 
all lattice sites for a particular realization of the disor- 
der (sample)u. A detailed numerical investigation reveals 
that the multifractal properties of the local Edwards- 
Anderson (EA) order parameter remains non-trivial far 
below the critical temperature. This contrasts with pre- 
vious investigations of the homogeneous pure ferromag- 
netic Ising model on the same lattice, for which the mul- 
tifractal spectra of the normalized local magnetization 



survives only at the critical temperatureo'B. Concerning 
the critical behavior of the Ising spin glass model, it was 
possible to obtain strong evidences of universal critical 
exponents and a unique value for the temperature of the 
critical point when several probability distributions for 
the quenched coupling constants are consideredQ'0. It is 
worth to mention that the values of the critical temper- 
ature and of the critical exponents associated to the or- 
der parameter and correlation length for Ising spin glass 
model on diamond hierarchical lattice with graph fractal 
dimension cLf — 3 are surprisingly close to those obtained 
by numerical simulations for the model defined on cubic 
latticesa. p 

In a recent publication, LuckEi proposed a heuristic cri- 
terion to account whether the geometric fluctuations in 
coupling constants and fields are relevant or irrelevant 
to alter (or not) the universality class of the aperiodic 
systems with respect to the corresponding homogeneous 
one. Luqkjs criterion plays the similar role of that due 
to HarrisEJ, that holds when quenched disorder is intro- 
duced in ferromagnetic systems. In the latter, the critical 
behavior is changed provided the critical exponent associ- 
ated with the divergence of the specific heat of the corre- 
sponding homogeneous system is positive. More recently, 
an exact derivation of an analog of. the Luck's criterion, 
was presented for both the IsingEl and the Pottai3E£l 
models, defined on the family of diamond hierarchical lat- 
tice filled with aperiodic interactions dictated by a gen- 
eralized Fibonacci sequence Using the transfer matrix 
(TM) approach, AndradeEJ obtained exact numerical re- 
sults for the thermodynamic and critical properties of the 
aperiodic ferromagnetic Ising model defined on the fam- 
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ily of diamond hierarchical lattices (hereafter DHL), cor- 
roborating the change of universality class of the model 
when relevant geometric fluctuations are considered ac- 
cordingly to the Luck's like criterion properly adapted to 
cope with hierarchical lattices. We notice that the exact 
scale invariance symmetry is the crucial property of the 
DHL's underlying the methods which link the thermody- 
namic functions of lattices with successive generations, 
as occurs for the above mentioned real renormalizatktn 
group scheme as well as for transfer matrices approachlifl. 

In this work we are mainly interested to investigate the 
role of geometrical fluctuations on the multifractal prop- 
erties of the local magnetization of ferromagnetic Ising 
model on the DHL's. We focus our study on two distinct 
models, both defined on DHL's with the same graph frac- 
tal dimension but with distinct topologies leading to rel- 
evant and irrelevant fluctuations. To deal with the multi- 
fractal properties of the local magnetization we apply the 
methodology developed to study the multifractal spectra 
of the order.uarameter of the homogeneous ferromagnetic 
Ising modelotiLj and the spin-glass Ising modelQ on the 
DHL's. As we will show, the caracter of the fluctuations 
influences the nature of the changes observed in the sin- 
gularity spectrum of the local magnetization. 

In section II, we describe the model Hamiltonian, 
present the sequences governing the aperiodic interac- 
tions and discuss the relevance of the geometric fluctua- 
tions for both lattices^-accordingly with the appropriated 
criterion presented inliil. We also review the main fea- 
tures of the flow diagram and the phase diagram of the 
aperiodic systems. Section III is devoted to review the 
method to obtain the local magnetization and to discuss 
the multifractal spectra for special realizations of the ge- 
ometric fluctuations. Finally we discuss and summarize 
our main conclusions in section IV. 



II. THE MODEL HAMILTONIAN AND THE 
APERIODIC INTERACTIONS 

A general hierarchical lattice is constructed according 
to a recursive rule where some bonds of basic unit are 
successively replaced by the basic unit itself. The gen- 
eral basic unit of the diamond hierarchical lattices (6, p)- 
DHL, is formed by two roots sites connected by a set 
of p parallel branches, each one containing a series of b 
connected bonds. A given generation is obtained by re- 
placing all bonds of the previous generation by the basic 
unit, as sketched in Figure [l]. For this family of hierar- 
chical lattices, the graph fractal dimension is given by 
dp = 1 + logp/log6, while its size, number of sites and 
the number of bonds are respectively given by L = £> , 
N s = (b - l)p[{bp) N - l]/(6p - 1) and N B = (bp) N , N 
being the number of generations. 

The reduced Hamiltonian for the system is 

- H/k B T = ^2 Ktjcncrj , (1) 

(i,3) 



where <jj = ±1 are the Ising variables and the sum 
runs over all pairs of nearest neighbors spins , 
Kij = Jij/ksT are the corresponding coupling constants, 
ks is the Boltzmann constant and T is the absolute tem- 
perature. 




(b-1 ) internal sites 



p branches 



FIG. 1. Basic unit of the general diamond hierarchical lat- 
tice (6,p)-DHL with p parallel branches each one with 6—1 
connected internal sites. 

In this paper, we consider geometric fluctuations in- 
troduced by the generalized Fibonacci inflation rule ex- 
pressed by 

(A, B) ^ (AB b -\ A b ) . (2) 

When applied to define the coupling constants of the 
hierarchical lattice (Jy = Ja\Jb), the above rule cor- 
responds to the following procedure: (a) a bond with 
coupling constant J a, ( J a > 0) should be replaced by 
a basic unit with p parallel branches each one contain- 
ing the first bond with coupling constant J a and the 
next (6 — 1) ones with coupling constants Jb; (b) a bond 
Jb,{Jb > 0) should be replaced by a basic unit with 
all bonds with coupling constants J a- This correspon- 
dence is displayed in Figures 2a and 2b, for the cases of 
the lattices with p = b = 2 and 3, respectively. After N 
steps, the length of the aperiodic sequences on which we 
base the present definition of the geometric fluctuations 
is given by b N , which corresponds to the number of bonds 
within each one of the shortest paths connecting the two 
root sites and measures the size of the lattice. We also 
notice that there are p N of such shortest paths, all of 
them with no common bonds and showing the same ape- 
riodic sequence. Therefore the systems can be regarded 
as formed by identical layers with the bonds J a and Jb 
arranged accordingly to the sequence (||). The number 
of distinct letters or coupling constants within two sub- 
sequent generations of the sequence are related by the 
substitution matrix M given by 
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whose eigenvalues are Ai = b and A 2 = (1 — 6). The wan- 
dering exponent which measures, the fluctuation in the 
distribution of elements within the sequence, is defined 
as 



log I A 2 
logAi 



(4) 



The sequences defined by (g) are said to be non-Pisot if 
oj > (b > 2), which means that the fluctuations are un- 
bounded, whereas the case u> — (6 = 2) corresponds to 
marginal case. On the other hand, for the Ising model on 
DHL with aperiodic interactions arranged by sequences 
such that Ai = 6, the geometric fluctuations are said to 
be relevant provided 



LO > LU C = 1 — 



(5) 



v being the critical exponent associated .with the corre- 
lation length of the homogeneous systemEJ. 

In the present work we consider two models: model 1 
is defined on a (2, 2)-DHL while model 2 is defined on a 
(3,3)-DHL. 





(a) 





(bi 



FIG. 2. Basic units of the (a) (2, 2)-DHL and (b) 
(3, 3)-DHL indicating the distribution of couplings J a and 
Jb accordingly with the generalized Fibonacci sequences. 

The corresponding lattices have been choosen to have 
the same fractal dimension (dp — 2), but the geomet- 



ric fluctuations are respectively irrelevant and relevant 
respectively accordingly, to eq.((|), as has been demon- 
strated by Pinho et altil. In the rest of this section we 
review the basic properties of the renormalization flow of 
the model Hamiltonian in these lattices. 

Consider a general (6,p)-DHL with the coupling con- 
stants defined by the inflation rule given by (^). Within 
the real space renormalization group scheme, the dec- 
imation procedure is carried on by partial tracing on 
the spins introduced in the last generation. This leads 
to exact scaling relations for the coupling constants J a 
and Jb, given respectively the following renormalization 
equations: 



t' A = tanh[ptanh 1 (t A t b I f 1 )} 
t' B = tanh[ptanh _1 (i^)] 



(6) 



where t x = tanh-Kx, x = A or B. 

The renormalization flow diagram in the (t^, tg) space 
has been carefully studied ir£lE3. The equations (g) have 
three fixed point solutions, all of them located along the 
manifold t A = ts'- = 0, = 1, and t*. Therefore, 
they are also the solutions of the corresponding homo- 
geneous system. The first two solutions refer to stable 
fixed points corresponding to the infinite (paramagnetic) 
and zero (ferromagnetic) temperature phases respectively 
while the third one is associated with the critical point 
governing the transition of the corresponding homoge- 
neous system. For the particular models considered here, 
these critical point solutions are exact and given by, 



t 



cl 



-(a c - 2/a c 



1) = 0.543689... for p = b = 2 (7) 



1) = 0.618033... for p = b = 3 



where a c — (3y33 + 17) 3. As shown i rfl the renor- 
malization flow space has two regions corresponding to 
the basins of attraction of the stable fixed points. The 
frontier of these regions, containing the critical point of 
the homogenous system, determines the phase diagram. 
However, as pointed out irclrO, linearization about the 
non-trivial fixed point indicates that it is a saddle-point 
for lattices with b — 2, (10 = < uj c ), while it be- 
comes fully unstable for lattices with b > 2, (uj > uo c ). 
Therefore, for model 1 the critical behavior should be 
governed by the fixed-point of the homogeneous system 
and the geometric fluctuations are irrelevant with re- 
spect to change the universality class of the model. On 
the other hand, for model 2, the non-trivial fixed point 
is fully inaccessible indicating that the geometric fluc- 
tuations are indeed relevant to change the critical be- 
havior. To grasp the critical behavior of the relevant 
aperiodic model one should explore further the flow di- 
agram gijzejn by eqs.(|^). This was already done by Had- 
dad et al!3 investigating the Potts model on the (3, 2)- 
DHL. Actually, whenever the fixed-point becomes fully 
unstable, they found the appearance of a two-cycle at- 
tractor associated with the novel critical behavior. For 
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the present model 2, this two-cycle attractor is given by 
coordination points (t* Al = 0.95116..., r t^ 1 = 0.46135...) 

0.99915. 



and (t* A2 = 0.54822..., t 



B2 



In Figure 3, we show the (tA, ts) flow diagram for the 
p = b = 3 case. Note that thetwo-cycle attractor is 
a saddle-point being unstable towards the trivial stable 
fixed-points and stable in the perpendicular direction. 
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A Two-cycle atractur 




FIG. 3. The flow diagram in the tA X f,s space for the 
real space renormalized couplings for model 2. The arrows 
indicates the flow of the second iterates. F and P denote the 
stable fixed-points corresponding to the ferromagnetic and 
paramagnetic phases. C labels the fully unstable fixed-point. 
A locates the two-cycle attractor and o the critical points 
obtained by the transfer matrix method, x labels the critical 
points obtained by inverting the renormalization equations. 
A and B give the frontiers of the basin of attraction within 
the manifolds ts — 1 and tA = 1 respectively 

The frontier between the basin of attraction of the sta- 
ble fixed-points contains both the two-cycle attractor and 
the fully unstable fixed-point. This frontier intersects the 
lines tg = 1 and tA = 1 in two points, respectively A 
and £?, and it manifold is associated with the critical tem- 
peratures of the model for different choices of the ratio 
R = Jul J A- 

The flow diagram can be numerically obtained by in- 
verting eqs.([3]) taking into account that the inverse equa- 
tions share the same trivial, non-trivial fixed points and 
two-cycle attractor as eqs.(|), but with inverted unsta- 
ble to stable manifolds and vice-versa. So, the iterates 
of these inverted equations starting close to the critical 
fixed point are now pushed towards line ACB, so that a 
relation ts = g(tA) can be numerically obtained ending 
to the phase diagram. 

In Figures 4a and 4b, we draw such relation in the 
(T/Ja, R) phase diagrams of both models, indicating the 
critical temperature obtained for some .special values of 
the ratio R by using the TM formalismEJ, as well as the 
corresponding exact values for the homogeneous mod- 
els (R = 1). The critical exponent v associated with 
the correlation length of the corresponding homogeneous 
systems can be exactly calculated as 
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FIG. 4. Phase diagram KbT/Ja X Jb/Ja for (a) model 1, 
(b) model 2. O (•) indicate numerical (exact) values obtained 
with the present methodology , while A indicate values ob- 
tained by transfer matrix method. A label the exact values 
for the homogeneous systems 



where 



dt'jt) _ Abptfjl + tfy^jl-tfy- 1 
~dT t=t , ~ ((l + t?)P + (l-tf)py 

and t'(t) is given by eqs.(||) when tA =ts =t. 



, (9) 



TABLE I. 


Exact values for v, 


ui and ui c for models 1 and 2. 




model 1 


model 2 


V 


1.338265788.. 


1.779416044... 


U! 





0.630929753... 


Ulc 


0.252764279.. 


0.438017880... 



In Table [jj we show the values of v, to and ui c for models 
1 and 2. It corroborates that for model 2 we have u> > 
lo c indicating that the geometric fluctuation should be 
relevant for this model. 
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III. THE LOCAL MAGNETIZATION 

The method of evaluation the local magnetization is 
based on the assumption that the model Hamiltonian for 
a lattice with N generations is equivalent to a reduced 
effective Hamiltonian of a single basic unit introduced in 
the N th step plus effective fields acting on the spins of its 
external sites and an effective interaction coupling these 
spins. This assumption, which has been proved to be for- 
mally correct for the homogeneous system, also holds for 
the present model, since no special condition is imposed 
to the coupling constants of the Hamiltonian. These un- 
known local effective interaction and fields represent the 
influence of the remaining lattice spin couplings trans- 
mitted by the external spins of that basic unit. The lo- 
cal magnetization of both internal and external sites of 
the basic unit of the effective system are calculated as 
a function of coupling constants {J a and Jb) and the 
unknown effective coupling and fields. Eliminating the 
unknown variables we end up to recursive relations for 
the local magnetization of the internal sites in terms of 
the corresponding values of the root sites. Now, consid- 
ering a N generation lattice and successively decimating 
the spins up to the first generation (basic unit) one can 
recursively calculate the local magnetization by choosing 
appropriated values for the magnetization of root sites, 
accordingly with the configuration of the phase of the cor- 
responding fixed point. InP the details of the calculations 
of the recursive equations-for the uniform ferromagnetic 
case is presented, while infcl these equations were general- 
ized for basic unit with arbitrary interactions in order to 
investigate the short range Ising spin glass model. The 
recursive equations for the latter can be straightforwardly 
applied to deterministic aperiodic Ising models express- 
ing the site magnetization of internal sites (a) in terms 
of the corresponding values for the external sites. The 
recursive expression for model 1 is given by : 



(a) = A\ (m) + B 1 (//) 



where 



Ax 
Si 



^a(1 — tAts) 

(1-4*1) ' 

*fl(l 



(10) 



(11) 



(1-4*1) 



for the (AB) non-homogeneous basic unit while 



A 1 =B 1 = 



tA 



(1 + 4) 



(12) 



for the (A, A)-homogeneous basic unit, as specified in 
Figure 2a. 

On the other hand, for model 2 we have two sets of 
recursive equations obtained for each one of the basic 
units shown in Figure 2b. Each set has two equations 
corresponding to the internal sites of each basic unit. For 



the non-homogeneous basic unit (ABB) these equations 
are given by 



(<7i) = A, (ji) + Bi (//) 



where 



An 



B 2 = 



t A (l 



L A L BJ 



(1 

4(1 ~ 4) 

(i-4*S)' 
t A t B (i - 1%) 
(i - 4* 



(i = l,2), 
(13) 



(14) 



4(1 



2 f 4 \ • 

2 f 2 \ 
A l BJ 



while for the homogeneous basic unit (AAA) they are 
expressed by 



{<r t ) = a <M> + A (//) 



where 



d=D 2 



D 1 = C 2 



t A {l-t A ) 



(1- 
4(1 



4) ' 

"4) 



(i = l,2), 

(15) 

(16) 



(1-4) 



To evaluate the correct value of the local magnetization 
m% = (ai) of each site, it is necessary to keep track of its 
local interaction configuration, i.e., to know with which 
choice of J a and Jb a given inner spin interacts with 
its root sites. In order to breakdown the global up-down 
symmetry and obtain non-zero values for the spontaneous 
magnetization, it is enough to apply a field on a given site 
of the lattice and let this local field equal zero at the end 
of calculation. However, within the present methodology 
this breakdown is easily done by assuming boundary con- 
ditions to root sites accordingly with the corresponding 
configuration of the stable fixed-point phase. This latter 
procedure is followed in the present work. 

Now we focus our attention to model 2 looking over 
the effects of the relevant geometric fluctuations with re- 
gard to the homogeneous pure model on the same lat- 
tice. First, we consider the local magnetization of the 
sites along one of the shortest path connecting the root 
sites. These values are assigned to points within the in- 
terval [0, 1], since the DHL is a fractal graph where the 
bonds mean " chemical distances" without any geometric 
sense. The resulting plot gives a representative picture 
of the local magnetization distribution, hereafter called 
magnetization profile, since all shortest path are equiva- 
lent as far the DHL graph symmetry is concerned. Figure 
5a shows the magnetization profile for the homogeneous 
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model 2 at the critical temperature of the correspond- 
ing saddle-point fixed-point, and Figures 5b and 5c show 
the profiles for the aperiodic model 2 at the critical tem- 
peratures of each one point of the two-cycle attractor 
for (which correspond to fluctuation parameter values 
Rd = Jb/Ja = 0.2706... and R c2 = J b /Ja = 6.3104.... 
For the first profile, in any of the renormalization steps, 
the coefficients ( 14[lq) have their values fixed by the cor- 
responding value at the critical point. For the other 
two profiles these coefficients alternate between the corre- 
sponding values at the two-cycle attractor. For values of 
J A and Jb other than those corresponding to the critical 
sets, the present methodology produces the local mag- 
netization as function of the temperature, but the above 
mentioned coefficients evolve under the renormalization 
process. In any case, the magnetization of the global root 
sites are held fixed with initial values equal to one, which 
corresponds to the ferromagnetic configuration. 

All profiles shown display many singularities and are, 
in some sense, self similar. However those of the aperiodic 
system clearly do not show the same symmetries as the 
one frem the homogeneous model (R = 1). As pointed 
out incru, for the homogeneous model the pattern sym- 
metries are solely related with the topologic symmetry of 
the lattice, imposed by the distribution of coordination 
numbers along the profile. However, for the aperiodic 
systems, the self similar symmetry of the lattice topology 
is superimposed by symmetries of the aperiodic sequence 
defining the distribution of interactions, so that the re- 
sulting pattern expresses the influence of both symme- 
tries. We emphasize that this feature does not depend 
on the relevant-irrelevant character of the fluctuations, 
and is observed also for the aperiodic situation of model 
1. It is worth to comment that for the spin-glass model 
the symmetry induced by the lattice topology is com- 
pletely washed out by .-the random quenched disorder of 
the coupling constants!! 

The profiles in Figures 5b and 5c are, at first sight, 
rather different, but a close analysis uncover the fact that 
they share the same structure. We note also the presence 
of two different scales corresponding to the regions where 
the letters A and B are more abundant. As the values 
of J a and Jb are rather distinct for the two points of 
the cycle, they enhance or depress the values of the local 
magnetization in these regions, leading to the two differ- 
ent shapes. 

It is important to observe that both profiles were ob- 
tained after eight generations, so that, in the flow dia- 
gram, the trajectories start and end at the same point of 
the cycle. 

If we consider an odd number of generations the form of 
the profiles would be reversed, indicating that the profile 
does not converge to a single form, but rather oscillates 
between two distinct patterns which reflect the properties 
of each of the points in the cycle. The same behavior is 
observed for every profile, evaluated with any choice of 
values for J a and Jb, at the corresponding value of T c . 
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FIG. 5. Magnetization profiles patterns of the model 2 
for lattices with TV = 10 hierarchies at the criticality. (a) 
R = 1 homogeneous system, (b) R c i = 0.2706... and (c) 
Rc2 = 6.3104... aperiodic systems at the fixed-points of the 
two-cycle attractor, respectively. 

As this alternating form of the profile pattern results 
from the different properties of the two points of the cy- 
cle, it is not observed in the case of irrelevant fluctua- 
tions. Indeed, for any value of J a and Jb, the trajectory 
at the corresponding T c is controlled by the properties of 
the critical point of the homogeneous system: the profile 
converges to a definite form, with a lower symmetry than 
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in the case of the homogeneous system. _ 
As previously observed for the homogeneous casesQ, 
the high degree of singularities suggests the multifractal 
analysis as a valuable tool to investigate the structure of 
the magnetization profiles for the models with R =/= 1. 
The F(a) functions of the measure defined by the nor- 



malized local magnetization i^i(L) 



were calculated by the method due to Chhabra and 
JensenliU which is based on a one parameter family of 
normalized measures, 



Mi(£) 5 



£>,-(£)]< 



(17) 



In Figures 6a and 6b we show the plots for both mod- 
els comparing the F(a) spectra of the aperiodic and the 
homogeneous system at criticality. Both curves have the 
same maximum value (F(a) max =l) since the measure 
has been arbitrarily assigned to the interval [0,1], that 
is, to the same fractal set (support) of dimension d = 1. 
Therefore, one must focus our attention to the values of 
the Holder exponent a. The minimum and the maxi- 
mum values of a reflect, respectively, how the measures 
of the most concentrated and the most rarefied inter- 
vals scale with the box width. For the homogeneous sys- 
tems the domain of the F(a) spectra can be analytically 
calculatecfil Moreover as was demonstrated iiillLj there 
is a straightforward relation between the Holder exponent 
and the critical exponents given by 



(18) 



where d is the fractal dimension of the support and (3 a 
is the critical exponent associated with the average mag- 
netization of the subset of sites whose the measure van- 
ishes with exponent a, as L — > oo. Since the most con- 
centrated measures of both models have finite values at 
the thermodynamic limit induced by the boundary condi- 



tions we have (3 ar 



and ou 



d — (3/v. Therefore 



oimin is solely related with the exponents describing the 
critical behavior of the whole system and should follows 
the rules for the universality class of the systems. More- 
over, for the subset of sites where the measure behaves 
with the Holder exponent a — d, the local magnetization 
behaves with the same critical exponent of the average 
magnetization of whole lattice, that is, f3^ = [3. This 
particular subset should also follow the rules of the uni- 
versality class, and therefore the point (d, F(d)) on the 
F(a) function should remains fixed when irrelevant geo- 
metric fluctuations are considered. 

From Figure 6b, it becomes clear that the aperiodic 
systems with relevant geometric fluctuations belong to a 
distinct class of universality regarding the homogeneous 
system, as the F(a)-function of the former is shifted to 
higher values of a. Moreover, we notice two distinct F(a) 
spectra, corresponding to distinct profiles with R = 1 and 
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FIG. 6. F(a) spectra of the normalized magnetization hier- 
archies at the criticality. (a) model 1 for lattices with N = 20 
and q € [—60, 60]. * labels the spectrum of the corresponding 
homogeneous systems, O the R = 2 aperiodic system and □ 
the R = 0.5 aperiodic system . (b) model 2 for lattices with 
N = 15 and q £ [—40, 50], * labels the spectrum of the corre- 
sponding homogeneous systems, Q an d D the spectrum of the 
R — 2 and R = 0.5 aperiodic system respectively . + and x 
label the spectrum of the R cl = 0.2706... and R c2 = 6.3104... 
fixed-points of the two-cycle attractor, respectively. 

On the other hand, as displayed in Figure 6a the F(a) 
functions of model 1 with R = 0.5 and 2 suffer minor 
changes due to finite size effects when compared with the 
one for R = 1. To handle with such differences we refine 
our calculations by doing a finite size scaling estimation 
of the values of a and F(a) for higher \q\, and by pro- 
ceeding a finite size scaling estimation of the exponent (f> 
governing the average magnetization at the critical point 
m ~ which by it turn is related with the critical 
exponents by <p = P/v- 

In Tables |0| and III, we present these estimations for 
both models respectively, comparing with exact values 
whenever is possible. We notice the slight difference 
when the a m i n of the aperiodic system of model 1 is 
compared with the exact value for the homogeneous sys- 
tem, indicating conservation of universality (see Table 
||). However, for model 2, qualitative differences are ob- 
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served when the a m i n values of the aperiodic and the 
homogeneous models are compared. In this latter case 
the values of a m i n oscillates as we go from N even to 
N odd lattices. As pointed out in our former discussion, 
these oscillations are a manifestation of the fact that the 
profile oscillates between two different patterns dictated 
by the properties of each of the points of the cycle. In 
Figure 7 we present the scaling for the a m i n of the F(a) 
function of model 2 for some values of the fluctuation pa- 
rameter R. For R c i and R c i these oscillations are clearly 
seen as we go from N even to N odd lattices. For fi- 
nite lattices and finite values of q, the oscilations on the 
scaling plot introduce significant errors bars for theval- 
ues of the corresponding exponents. However, as N goes 
to infinity, we expect that a m in converges to a unique 
universal limit for any value of the fluctuation parameter 
other than R = 1. 

ob- 



In Table III, we also compare the values of a r 



tained within the present approach with the ones deduced 
from the results for (3/u, calculated using the transfer 
matrix approach developed by one of ualj. It is worth 
to mention that this latter approach allows to work with 
higher values of N, ending up with more accurate values 
for the exponents. 



A 
A g 
A $ 
A ^ 



A R=1 

+ R c „=0.2706. 

X R — ,=0.6310. 

□ R=0.5 

O R=2 



N In 3 



FIG. 7. Scaling of the a m in of the _F(a)-function of model 
2 for values of R = 1, R = 0.5, R = 2, R cl = 0.2706... and 
Rc2 = 6.3104.... The corresponding exponents are given on 
table III. 



IV. SUMMARY AND CONCLUSIONS 

We investigate the role of the geometric fluctuations 
on the local magnetization and critical properties of the 
ferromagnetic Ising model when the coupling constants 
are defined by deterministic aperiodic sequences. Two 
distinct models were considered as a matter of compara- 
tion. Model 1 (model 2) was build with irrelevant (rele- 
vant) Fibonacci like geometric fluctuations with respect 
to changes in the critical properties. The models are de- 
fined on different diamond hierarchical lattices but with 
then same fractal dimension dp = 2. The phase dia- 



gram of both models were obtained and carefully ana- 
lyzed. The local magnetization of both models was cal- 
culated by an exact recurrence procedure as function of 
the temperature but we concentrate our analysis close to 
the critical point. When the relevant (irrelevant) geo- 
metric fluctuation is concerned the pattern of the local 
magnetization profile exhibit a strong (weak) change of it 
symmetry as compared with the one of the corresponding 
homogeneous system. For the model 1, the self-similar 
pattern induced by the lattice topology is predominant, 
as it occurs for the homogeneous system. 

However, for model 2, the pattern symmetry of the ho- 
mogeneous system is strongly affected by the one dictated 
by aperiodic distribution of coupling constants and the 
local magnetization profile resulting with a quite different 
distribution of singularities. This break of symmetry is 
straightforwardly related with the change in critical prop- 
erties, as revealed by the multifractal analysis of the mag- 
netization profiles at the criticality. The F'(a)-functions 
of the model 1 is quite close to the one of the correspond- 
ing homogeneous systems for all values of domain of the 
spectrum, reflecting the same critical behavior of the lat- 
ter. However, the domain of the i ? (a)-functions of the 
model 2 is shifted to higher values of a, with respect to 
the one of the corresponding homogeneous systems, lead- 
ing to change the universality class of the model regard- 
ing the relative values of the ferromagnetic coupling con- 
stants. The critical exponent associated with the magne- 
tization of model 2 were calculated by means of the F(a)- 
function as well as by scaling the average magnetization. 
The results are in good agreement with the ones obtained 
through the thermodynamic functions calculated using a 
quite different approach previously developed by one of 
us. Therefore, we have demonstrated, by means of exact 
numeric calculations of the local magnetization of useful 
models, that as far as relevant deterministic geometric 
fluctuations is concerned the resulting disorder in the dis- 
tribution of aperiodic ferromagnetic interactions leads to 
a strong break of symmetry in the patterns of the local 
magnetization following the change of universality class. 
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TABLE II. Domain of the F(a) spectra for the aperiodic and homogeneous systems for model 1. The numerical values 
obtained by the scaling of a m in and a max of the aperiodic systems were calculated for lattices with N=10 to 22 and for 
q — ±50. The finite size scaling of the average magnetization was performed for lattices with L = 2^, N — 10 to 22. 



method 




R=l 


R=2 


R=0.5 


analytic a 




0.8791464205 






scaling b 




0.87895(6E-5) 


0.8793(9E-4) 


0.876(E-3) 


scaling 6 


l-P/v 


0.87900(3E-5) 


0.8794(9E-4) 


0.879(E-3) 


transfer matrix c 


1-P/v 


0.87915(E-5) 






analytic a 




1.04604781... 






scaling b 




1.03920(7E-5) 


1.0467(6E-4) 


1.0474(5E-4) 



a RcJ3Q 



b pres«iit work 
c Ref.y 



TABLE III. Domain of the F(a) spectra for the aperiodic and homogeneous systems for model 2. The numerical values 
obtained by the scaling of a m in and a m ax of the aperiodic systems were calculated for lattices with N=6 to 15 and for q = ±40. 
The finite size scaling of the average magnetization was performed for lattices with L — 3 N ,N = 6 to 15. 



method 




R=l 


R=6.3104... 


R=0.2706... 


R=2.0 


R=0.5 


analytic a 




0.8760357586... 










scaling b 




0.8768(5E-4) 


0.95(E-2) 


0.94(E-2) 


0.9314(7E-4) 


0.9219(8E-4) 


scaling b 


1-P/v 


0.8768(5E-4) 


0.95(E-2) 


0.94(E-2) 


0.9301(6E-4) 


0.9178(6E-4) 


transfer matrix c 


1-/3/!/ 


0.8760(E-4) 


0.95481(E-5) 


0.955(E-3) 






analytic a 




1.068947633.. 










scaling b 


Otmax 


1.0643(6E-4) 


1.1660(5E-4) 


1.162(5E-3) 


1.166(5E-3) 


1.160(5E-3) 



" Rcffl 



b preseint work 
c Ref.y 
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